Disordered and Ordered States in a Frustrated Anisotropic Heisenberg Hamiltonian 
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We use a recently proposed perturbative numerical renormalization group algorithm to investigate ground- 
state properties of a frustrated three dimensional Heisenberg model on an anisotropic lattice. We analyze the 
ground state energy, the finite size spin gap and the static magnetic structure factor. We find in two dimensions a 
frustration-induced gapless spin liquid state which separates two magnetically ordered phases. In the spin liquid 
state, the magnetic structure factor shows evidence that this state is made of nearly disconnected chains. This 
spin liquid state is unstable against unfrustrated interplane couplings. 
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Low dimensional quaritum magnets are currently the object 
of an important interest 1 1]. This stems from the rich physics 
which is displayed by these sytems due to their reduced di- 
mension and competing interactions which often push the 
transition to ordered states to very low temperatures or even 
preclude the onset of magnetism at all. Geometric frustration 
is one of the effects which are believed to lead to possible non- 
classical states. A non-classical ground state exists in a pure 
one-dimensional quantum antifeiTomagnet. The ground state 
is disordered and the low energy excitations are spinons which 
carry fractional quantum numbers. A fundamental question is 
whether spinons can survive in higher dimension. 

The relevance of these questions has been substantiated in a 
recent neutron scattering experiment |2] reported in the quan- 
tum magnet CS2CUCI4. This compound is a quasi-two dimen- 
sional spin one-half Heisenberg antiferromagnet. It is made 
of anisotropic triangular planes (with J' = J" « 0.34J in 
EqQ) which are weakly-coupled by an exchange which is 
roughly — 0.045 J. An incommensurate Neel state is sta- 
ble below Ttv — 0.62if at zero magnetic field. In the two 
dimensional regime above T/v, the dynamic correlation dis- 
plays a highly dispersive continuum of excited states which is 
a signature of spinons. This finding together with the fact that 
J' is in the same order as J led to the conclusion that above 
Tn, CS2CUCI4 is a 2D spin liquid with fractional excitations. 

In order to study these effects of frustration and dimen- 
sional crossover, we will consider the following anisotropic 
3D Heisenberg Hamiltonian: 

i,l^k i,l 
J" Sij^k^i+i^i+i^k + Jz ^ ^i,l,k^i,l,k+li (1) 

where the i index represents sites, the I index chains and 
the k index planes, the exchange couplings are such that 
< J = 1 and J', J" < 1. When J' ^ and J" ^ the 
system is frustrated. In this case the Quantum Monte Carlo 
method, which is so far the most reliable method of investi- 
gation of spin systems, is plagued by the sign problem. Al- 
ternative approaches |3, 4] have investigated the possibility of 
a spin liquid state in the 2D regime ( = 0) of the Hamil- 
tonianQ. In Ref. |3], this 2D model was studied using a 



combination of Ising and dimer series expansions. Consid- 
ering the extrapolation of their results to the limit of weakly 
coupled chains, the authors was not able to conclude clearly 
due to the limitation of their technique. They argued that ei- 
ther a spiral ordered or a nearly critical disordered phase are 
possible in this regime. The dynamical susceptibility of this 
model was computed within the Random Phase Approxima- 
tion (RPA) |4| using an essentially exact expression for the 
ID chain dynamical susceptibility. For the 2D system incom- 
mensurate order with exponentially small characteristic wave 
vector is predicted. But this prediction was inspired by the ex- 
perimental observation. It is impossible to tell from this study 
if this is an intrinsic behavior of the HamiltonianlO). 

We wish to present in this letter an ab-initio computa- 
tion of the ground state static properties of the Hamilto- 
nian([0. For this purpose, we use a recently proposed per- 
turbative density matrix renormalization group (DMRG) ap- 
proach This perturbative DMRG method which 
has been so far used for only anisotropic 2D systems is ex- 
tended here to anisotropic 3D systems. We will explore the 
following issues: (i) what is the nature of the ground state in 
the 2D regime of the Hamiltonian Q? (ii) Is a spin liquid 
state favored when J' = J" ~ 0.3 as found experimentally? 
If so, what is the nature of this spin liquid state? (iii) Does Jz 
restore a magnetic state starting from this eventual spin liquid 
state? If so, is this order incommensurate? 

The perturbative DMRG which will be used in this study is 
a particular case of a more general matrix perturbation method 
based on Kato-BIoch expansion 1^ fioll which was recently 
introduced by one of us 1 7]. In the first step, the usual ID 
DMRG method 1,1 1,1 is applied to find a set of low lying eigen- 
values e„ and eigenfunctions |0„) of a single chain. One 
can note that during this step, if one wishes to study not too 
large lattices, it is preferable to use the exact diagonalization 
method instead of the DMRG. The advantage of the exact di- 
agonalization method is that by selecting both the total spin S 
and the momentum in addition to the Sz component of S as 
done in DMRG, it will lead to a better estimation of the low 
energy Hamiltonian of a single chain. 

In the second step, the 2D Hamiltonian (Eq.Q with = 
)is then projected onto the basis constructed from the tensor 
product of the |(/)n)'s. This projection yields an effective one- 
dimensional Hamiltonian for a single plane (we drop tempo- 
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FIG. 1: Static spin structure factor 5(fci:, fcj,) for La; X Ly = 24x25, 
J' = 0.3, J" = (upper graph) and J' = 0.3, J" = 0.6 (lower 
graph) 

raly the plane index k), 

["] ' 

J"^S,S,+i (2) 

I 

where is the sum of eigenvalues of the different chains, 
= Tli^ni^ |^[ri]) is the con-esponding eigenstate, 
l*[n]) = l</'ni)|<^n2)---l^ni,); the composite chain-spin op- 
erators on the chain I are S; = (Si;, S2/, ...S^;), where the 
running index labels sites in a chain I. Note that the products 
of spin-chain operators in Eq.|2lare different for the terms in- 
volving J' and J" . The term with J' couples spins with same 
site index on neighboring chains, when the one with J" cou- 
ples spins with site index i on chain I and i — 1 on chain I + 1 
are coupled. These renormalized matrix elements on the sin- 
gle chain basis are 



The effective Hamiltonian (|3 is one-dimensional and it is 
also studied by the DMRG method. The only difference with a 
normal ID situation is that the local operators are now mxm 
matrices, where m is the number of states kept to describe the 
single chain. 

In the third step, the same procedure is repeated in order 
to go from 2D to 3D once is set on. One then obtains the 
following effective ID Hamiltonian for the 3D system: 

^3 « E ^[P] l*bl ) ^* W \ + (4) 

[p] k 

where E'^^^ are the sum of eigenvalues of the decoupled 
planes, |^'[p]) are the corresponding eigenstate, and the 
composite plane-spin operators on a plane k are S'^ — 
{Sik, S2k, •■•Spfc), the running index labels chains in a plane 
k. The generation of the effective Hamiltonian is iden- 
tical to that of the effective Hamiltonian (|2ji described in 
Refs.|^0]. We target the spin sectors Sz — 0, ±1, and 
±2 to generate a low ener^ Hamiltonian which describes 
the 2D systems. In Refl*^, 8], this perturbative DMRG was 
tested against the stochastic series expansion quantum Monte 
Carlo (QMC) and exact diagonalization methods respectively. 
The agreement was very good for small transverse couplings 
and not too large lattices. These comparisons showed that the 
method is well controlled and its accuracy can systematically 
be improved by increasing the number of states kept during 
each step. 

We first set — Q and study the ground state properties 
of a single plane. We studied lattice sizes ranging from 8x9 
to 32 X 33 for J' = 0.15 and J" varying from to 0.3; then 
J' = 0.3 and J" going from to 0.6. We tipically keep nix ~ 
128 states in the single chain calculations; i.e. the whole chain 
is described by 4 x 128^ states. Among these states we keep 
only up to niy = 64 states in the second step; i.e. a plane is 
described by 64"^ states. The truncation errors are smaller than 
10~®. The reason for this abnormally small truncation error, 
which is related to the use of three blocks instead of four, has 
been discussed in Ref.|51. 

We found a similar qualitative behavior for these two sets of 
parameters. When J" < J', the correlation in the transverse 
direction are AFM; the static structure factor S{kx, ky) shown 
in the upper graph of Fig.fOfor a 24 x 25 lattice, displays a 
maximum at (kx, ky) = {it, tt) (Note that the triangular lattice 
is equivalent to a square lattice with a coupling along one diag- 
onal only. We thus adopt a square lattice notation throughout 
this study). We can thus conclude that in this regime, a Neel 
state is stable. Now if J" > J', the maximum in S{kx,ky) 
shifts to {tt, 0) as seen in the lower graph of Fig.^ This is a 
signature of the coUinear phase as in Ref.|6]. This behavior 
was expected and has been found in other frustrated models: 
if one of the competing exchange parameters is dominant, the 
coiTesponding ordered state ( the order which might prevail 
in the absence of the other competing exchange) is favored 
OBl- A disordered state can only be expected in the re- 
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a single chain (diamonds). In the inset, error in the binding energy 

5Eb = \Eb - Ei°\ (in log scale) as a function of J" for J' = 0.3. ^iq 3. ^(fc^^ fc^) fo]. x Lj, = 24 x 25, J' = 0.3 and J" = 0.3. 



gion of parameter space where the competing exchanges are 
close. 

We now concentrate on to the regime J" = J'. The pos- 
sibility of a disordered critical state was raised by the series 
expansion study of Ref|3]. The expansion becomes inacu- 
rate for small values of ,/'/ J, these authors could not reach 
this critical regime. Our data show the existence of such a 
disordered state for J' = J" = 0.15, 0.3, and 0.45. Fig. 13 
shows for J" — J' ~ 0.3 the finite size gap A2D of the 2D 
system and that of a single chain Aid. A2D is found to be 
always smaller than Au^. Since A1/3 when L ^ 00 
and m 00, we can thus conclude that the 2D system is also 
gapless in the thermodynamic limit. The ID version of this 
model, the two-leg zigzag ladder 1 14], has an exponentially 
small gap {Aid oc exp(— aJ'YJ) and a is a constant) hard 
to extrapolate from a finite size analysis. That small gap co- 
exists with ferromagnetic correlations between adjacent spins 
on the weakest bonds of the dimerized states. None of these 
concurrent signatures of the ID frustrated system are observed 
in the 2D model where the ferromagnetic correlations favor a 
collinear state. 

For J' = J", S{kx, ky) shown in Fig.|3lis structureless in 
the ky direction and it retains its maximum at tt in the k^ di- 
rection. This indicates that the chains are disconnected. The 
small bump which can be observed in S{kx: ky) near (tt, tt) is 
a consequence of very small short-range transverse AFM cor- 
relations. Thus the spin liquid state is mostly ID. Our conclu- 
sion of a disconnected chains ground state is also supported by 
the transverse local bond strength | (S^; Su+i) \ and the binding 
energy of the chains Eb = i?iD — E2d/{L + 1) (not shown), 
where Eid is the ground state energy of a chain of length 
L and E2D the ground state energy of a L x L + 1 lattice. 
|(S,;S,;+i)| decays from 0.0818 in the Neel state (J' = 0.3, 
J" = 0) and 0.0766 in the collinear state (J' = 0.3, J" = 
0.6) to 0.0132 in the disordered state (J' = 0.3, J" = 0.3). 

Indeed, one expects the perturbative DMRG method to be 
valid only at small couplings. One may thus question its valid- 



ity for relatively large couplings such as J' — J" — 0.3, and 

0. 45. In ReffSi] we have shown that, when larger values of the 
transverse coupling are used, the two-step DMRG is indeed 
less accurate when a genuine two-dimensional state emerges. 
But at the maximally frustrated point (which is at J' = J" ) 
the method remains surprisingly accurate even for intermedi- 
ate values of J' and J" . This is seen in the inset of Fig|2]which 
shows that 6Eb, the difference in binding energy between the 
perturbative DMRG and exact diagonalization, is minimum at 
J' — J" . This is because at the maximally frustrated point, 
the interchain correlations are very small. Thus consistent 
with the observation made on S{kx, ky), in the ground state, 
the eigenstates of an isolated chains are a good approximation 
for the 2D system. 

We now study the stability of the 2D spin liquid state 
against interplane perturbations. The 3D results are less ac- 
curate than the 2D ones. The maximum truncation error was 
about 10"'^ during the generations of the effective 2D Hamil- 
tonian to be used as the starting point in the 3D computations. 
This relatively large truncation error is due to the fact that five 
spin sectors are targeted. We thus were restricted to smaller 
2D lattices. The largest lattice we studied was 16 x 16 x 17. 
We set J' = J" = 0.3 and we vary J, from 0.05 to 0.15. 
Fig. |4] displays S{kx,ky) for = 0.2 in the middle plane, 

1. e. the 9*'* plane in a 16 x 16 x 17 lattice. One can observe 
the emergence of a small peak at (tt, tt) which is due to the in- 
terplane coupling. We find that as one may expect, the height 
of this peak increases with increasing J^. Since the coupling 
in the third direction is not frustrated, magnetic energy can be 
gained by coupling the planes. As a consequence, the spin liq- 
uid state which results from the inability of the chains within 
the planes to couple effectively is destroyed. 

To summarize, we have studied the ground state properties 
of a 3D anisotropic Heisenberg model. We can now address 
the issues mentioned in the introduction which were motivated 
by the experimental results of Ref.|2]. (i) In the 2D regime 
of this model, the magnetic structure factor, the binding en- 
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FIG. 4: S{k^, ky) for x Ly ^ IQ x 16, J' = 0.3, J" = 0.3 and 
= 0.15. 

ergy display the properties of a spin liquid state. This spin 
liquid state separate two regimes with long-range order, (ii) 
This spin liquid state is gapless and is made of nearly inde- 
pendent chains. It is remarkable that it survives even when 
the interchain exchange are of the same magnitude as the in- 
trachain exchange. Since it has been previously reported in 
the crossed chain model 1 15] and on a square lattice, it seems 
to be generic of models of chains coupled with a frustrating 
interaction, (iii) Turning on the interplane coupling restores 
magnetism. But we do not find any tendency to incommen- 
surate ordering. The discrepancy on this point between our 



simulation and the experiment f?) may be due to the fact that 
in the compound CS2CUCI4, consecutives triangular layers are 
sligthly rotated from each other. This effect can be taken into 
account by adding a Dzyaloshinskii-Moriya interaction into 
the HamiltonianQ. This term could be the source of the in- 
commensurate ordering. 

We thank Jim Allen for a critical reading of our manuscript. 
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